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Abstract 

We consider the possible production of stable lightest hrst level KK par¬ 
ticle (LKP) in baryonic gamma ray bursts (GRB) out flows. We numerically 
computed the energy-dependent cross-sections of Kaluza-Klein (KK) excita¬ 
tions for the Standard Model gauge bosons, 7 and Z”. Next, we determined 
the feasibility of producing these KK excitations in gamma-ray emitting re¬ 
gions of GRBs. We found that a GRB fireball that accelerates baryons to 
energies greater than 10’^^ eV could produce KK excitations out to approx¬ 
imately 10’^^ cm from the central engine, indicating that GRBs may be a 
signihcant source of the LKP. Finally, we explore the potential observational 
consequences of our results. 
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1. Introduction 

Kaluza-Klein (KK) theory allows all standard model (SM) fields to prop¬ 
agate into compact extra dimensions [T] , which is a simple extra-dimensional 
extension of the SM with a dark matter candidate. All helds are assumed 
to propagate in a flat space-time M4 x Si where is the extra dimension 
compactihed on a circle. Since the helds are conhned and must be continu¬ 
ous, a discrete series of particles arises. The momentum of the particles in 
the extra dimension appears as part of the apparent 4D mass. This creates 
a tower of KK resonances (or excitations) for each particle. Each excitation 
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is heavier than the last and has the same qnantnm nnmbers and spin as the 
4D particle. 

The mass of the nth mode in the KK tower, m”, is given by, 



where mP is the mass of the 4D particle (also called the zeroth order KK 
mode) and R is the radius of the compactihed extra dimension. 

The momentum conservation in the extra dimensions leads to the con¬ 
servation of KK number at each vertex. Conservation of KK parity requires 
that interactions only occur between even-number or odd-number KK modes. 
Thus, there exists a lightest hrst level KK particle (LKP), usually a linear 
combination the KK photon and the KK with suppressed decay cross sec¬ 
tion to the SM zero modes [2]. The unstable KK modes will quickly decay 
into the SM zero mode or the stable LKP. The identity of the LKP crucially 
depends on the mass spectrum of the hrst KK mode. Moreover, the LKP 
is considered a promising dark matter candidate because it is long-lived and 
has low annihilation cross sections [3]. The current Large Hadron Collider 
(LHC) limit on the compactihcation scale is R~^ > 715 GeV and the mass 
of the LKP Mkk > 1-4 TeV |1]. 

Experimentally, direct detection searches are attempt to isolate rare in¬ 
teractions between dark matter (DM) and ordinary matter. For example, 
the Large Underground Xenon (LUX) [5] and XENON-100 |6] dark matter 
experiments searched for events in liquid xenon target by recreating the light 
received by photomultipliers. These efforts are limited by a DM candidate’s 
mass and its scattering cross-section, and their maximum possible sensitivity 
depends on background signal from coherent neutrino scattering. 

There are also indirect methods that search for evidence of rare annihila¬ 
tion or scattering events elsewhere in the universe. In fact, the experimental 
results from numerous instruments, e.g. Fermi and PAMELA, have noted an 
excess in the cosmic electron-positron spectrum, which may be consistent 
with the annihilation events of 3 TeV DM [aiH]. However, these Endings 
are open to interpretation, especially since the signals cannot be localized to 
known DM sources |9]. 

On the theoretical front, recent work on indirect astrophysical detec¬ 
tion of exotic dark matter candidates ra explored radiation signatures from 
electron-dark matter scattering near active galactic nuclei, reaching the rather 
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pessimistic conclusion that such signals are probably undetectable with ex¬ 
isting and near future instruments. In the more specihc context of active 
galactic nuclei (AGN) jets, [TT] performed a more detailed calculation that 
included specihc candidate particles including the KK. They found that res¬ 
onances dictated by exotic particle mass spectrum may dominate scattering 
between proton (or electron) and dark matter candidates, resulting in a po¬ 
tentially detectable break in the gamma ray spectra from AGNs. However, 
these studies assumed that the exotic dark matter already exist in sufficient 
quantities, and did not consider how such particles may have been produced. 

In this work, we propose a potential astrophysical scenario for the produc¬ 
tion of the LKP as DM, and discuss the possible observational implication. 
DM searches that rely on rare interactions, restricts the number of events 
that we can observe. Particle accelerators overcome this by creating events, 
yet are limited by the energies that are required to create these events. How¬ 
ever, some cosmic particle accelerators in our universe that do not have the 
same energy constraints, including active galactic nuclei, supernova remnants 
(SNR), and gamma ray bursts (GRB). The improving sensitivity of gamma- 
ray instruments may enable us to search for potential DM matter signatures 
in the gamma ray regime from such cosmic accelerators. 

We show that a combination of the current theories for KK DM and GRB 
dynamics means that GRBs could produce KK DM in appreciable quanti¬ 
ties. GRBs are attractive possibilities because they are, at least in principle, 
plausible sources for ultra-high energy cosmic rays at energies between 
to eV, especially if the lower energy cosmic rays can come from other 
systems such as supernova remnants da [13]. In this context, baryon-loaded 
GRB hreballs may accelerate protons and neutrons to above 10^° eV (Hj, pos¬ 
sibly via a combination of hrst and second order Fermi acceleration [131 at 
shocks. Furthermore, GRBs create mildly relativistic shocks with the inter¬ 
stellar medium (ISM) [16]. They may also create strongly shocked regions, 
when hreball outflows, with randomly distributed energies, collide. More¬ 
over, the millisecond variability of GRB light curves may rise from multiple 
shells ejected by the central engine, so it is possible for there to be several 
internal collisions. Particle encounters about these shocks generate particle 
showers which are translated into the gamma-rays. These particle showers 
may be frequent and energetic enough to produce KK excitations. Beyond 
high energy cosmic rays physics, these hadronic models also have garnered 
signihcant attention due to their potential to explain a range of observations 
[niiiHiiin]. We refer interested readers to review articles such as |2ni[2ll[22] 


3 


and [23] as well as further references for more details regarding the rich and 
fascinating phenomenology of GRBs. 

We use the hreball model from m to estimate physical properties of the 
hreball. We assume a baryon-dominated hreball as such outhows are more 
likely to generate signihcant quantities of high-energy particle encounters 
compared to Poynting dominated scenarios. Next, we numerically compute 
the energy-averaged cross-section of processes involving KK particles, and 
combine these results with the collision region parameters to obtain an exci¬ 
tation rate and optical depth. These results then enable us to identify, at a 
qualitative level, potential signatures of the KK theory in GRBs. 

Our paper is organized as follows. In Section 2.1 we outline our numer¬ 
ical cross-section computations and some necessary aspects of KK particle 
physics. Section 2^ then explains the GRB dynamics relevant to our calcu¬ 
lation while summarizing relevant results from |23|. Section 2.3 the explains 
the reaction rate calculations. In Section we present our main hndings. In 
Section]^ we discuss our results and the potential observational implications. 


2. Methods 

2.1. KK excitation at high energies 

We used the computer program Pythia 8 [25| to simulate proton-proton 
collisions with and without a TeV“^ sized extra dimension [2]. We assume 
that any resulting KK tower quickly decays into either the LKP or SM par¬ 
ticles. The LKP was given a set mass of 1.5 TeV, slightly above the current 
limits set by the LHG. These calculations are consistent with ongoing dilep¬ 
ton searches at the LHG [3]. To represent all the relevant interactions, we 
included the processes // —j- /^*)n —t FF where f{F) can be any 

initial (hnal) state fermion. 

The KK bosons mediate processes between fermions [2S|- The KK ex¬ 
citation masses are dependent on the extra dimension size R through the 
relation m* = where m* is the mass of the KK excitation. Thus, for 
the LKP KK excitation masses and are given by, 

+ {n ■ m*Y (2) 

n ■ m* (3) 

We performed our calculations at logarithmically spaced energies starting 
from the LHG energy, 14 TeV, to the maximum astrophysical energies, ~ 10® 


(n) 

mz* = 

(n) 

ml, = 


4 






TeV. We report the computed cross sections cr{E) in figure where E is 
the beam energy. A larger cross section indicates a higher probability that 
fermionic processes will be mediated by a KK boson. All simulated values 
were for a 95% confidence interval. 


2.2. GRB Dynamics 

The excitation rate per unit volume is a Lorentz-invariant quantity and 
is dependent on the cross section for KK excitation and the number and 
energy of the protons in the outflow. In GRBs, for these KK excitations to 
be appreciable, they must proceed quickly with respect to the time-scales of 
the bursts. The gamma-ray flux from GRBs may last anywhere from two 
to several hundreds of seconds, with a few notable exceptions [27]. In the 
previous section, we calculated the probability that a collision will give a KK 
excitation. In this section, we describe the energy distribution of protons in 
the outflow. 

In our model, the central engine is spherical and its mass is concentrated 
in baryons. Since GRBs preferentially occur in areas with low metallicity 
|28j ■ we assume a flow composed of free neutrons and protons. For long 
GRBs, the progenitor object is believed to be a massive star collapsing into 
a blackhole, so that the initial size is on the order of the Schwarschild radius 
Rq 2GMq/c^ ~ 10^ — 10^ m. The energy of the burst is related to its 
mass by Mq = Eq/Tc^, where isotropic energy output of a GRB can reach 
up to 10^^ ergs [29]. The energy imparted to this mass Mq Eq/c^^ within 
the initial radius Rq, forces a relativistic expansion of the material. The 
relativistic outflows of GRBs have a coasting Lorentz factor on the order 
of To 100 — 1000. An average coasting Lorentz factor of 100 and an 
isotropic energy output of 10®^ ergs would suggest an outflow mass with 
Mo 10 ^Mq. The coasting Lorentz factor is reached at a saturation radius 
Rg. The Lorentz factor as a function of distance from the central engine R 
is therefore 


R/ Rq if R < Rsat 

To if i? > Rgat 


(4) 


where saturation radius is defined by Rg — toTo. 

As the fireball expands into a relativistic outflow, the initial thermal 
energy of the fireball is progressively converted to kinetic energy. From the 
observers frame, the outflow organizes into a radial shell of width Ai? ~ Rq 
with the mass of the outflow given by M = Nmi,, where N is the number of 
baryons and rub is the baryon mass. We assume that protons and neutron 
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have approximately the same mass, hence ~ rrip. The comoving proton 
density is given by, 




EqX 


4:/3nRQC^mbT 


(5) 


where x is the proton fraction, i.e. x = np/ribi which we assume to be 0.5. 

As the relativistic outflow expands into the ISM, it produces shocks which 
inject baryons from the ISM into the shell. GRBs generally occur in dense 
star forming regions, therefore we assume that the ISM has a number density 
of nisM = 1 cm“^. The shell will continue to propagate with Lorentz factor 
To in the coasting phase until it reaches the spreading phase at rs- In the 
spreading phase, the shell expands with its local sound speed due to the 
velocity difference within the shell. The shell continues to propagate and 
expand, reaching the deceleration radius r^, it begins to lose kinetic energy 
to the cold ISM. At which point. 


td 


3E 


0 




( 6 ) 


Internal collisions occur when a shell catches up to another shell. Consider 
a system where each of the shells ejected from the central object have random 
mass-energy injections and dimensionless velocity (3, where T = 1/v^l — (3'^. 
This phenomenon produces a number of internal collisions between these 
shells. 

Consider a faster shell ejected sometime after a slower shell, AT. The 
faster shell, henceforth denoted by the subscript /, will sweep into the slower 
shell, henceforth denoted by the subscript s. The faster shell injects energy 
into the slower shell, thus it will be called the injective shell. The slower 
shell receives the injection of energy in the form of baryons. The collision 
creates a shock between the two shells, similar to the shock created between 
the shell and the ISM. The radius at which these two shells collide is 


r„„ ~ 2T]c(^T)|(l - r./r,), (7) 

where T f and T^ are the Lorentz factors of the fast and slow shell respectively. 
Notice that in the limiting case where Tj 3> T^ that this reduces to Vcoii — 
2r^(cAT). For a shells with Tj = 300 and T^ = 1, rcou ~ 10® m if the fast 
shell is ejected 10 ms after the slow shell. 
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2.3. Excitation Rate and Optical Depth 

KK excitation may occur in shocked regions as the outflow collides with 
either the ISM or a slower shell. The shocks facilitate diffusive and stochastic 
acceleration of the baryons, which results in a power-law energy distribution 
np{E) oc {E/Eo)~°', where the spectral index is a ~ 2. We take the width of 
the shocked region to be AR ~ Rq. 

Once the protons have been accelerated above the resonant energy, col¬ 
lisions may produce KK excitations. To quantify the importance of KK- 
mediated processes, we compute an excitation optical depth 


tkk 


nb{r){aKK{r))AR 


( 8 ) 


as a function the distance r where the shell smashes into either a slower shell 
or the ISM. Here we take the region where significant KK mediated collisions 
occur to be about the shell width AR, and Ub is the baryon number density 
in this region that depends on both r and the nature of the medium that the 
shell runs into. To compute tkk, we use a energy-averaged cross section 


{(^KK{r)) = 


Ja{E,r)(p{E,r)dE 
I4>{E,r)dE ’ 


(9) 


where (j){E,r) = nb{E,r)v ~ nb{E,r)c is the energy-dependent proton flux 
in the highly relativistic flow. Here we assumed that nb{E,r) is approxi¬ 
mately spatially constant across AR. We report tkk in figure for different 
scenarios. 

Finally, we calculate the production rate. The nuclear collision rate is 
given by, hcoii{E) = np^i{E)a{E)np^ 2 clT The protons in 77,^2 are as¬ 
sumed to be stationary with respect to np^i{E). This assumption is intu¬ 
itive for shocks with the ISM. While injections of protons by the ISM may 
be accelerated in these shocks, it is much more reasonable that they collide 
with the already accelerated protons. Thus, the production rate is nKK{E)- 
= np^s{E)a{E)npjsMc/Trei, where Up^g and UpjsM are the particle densities 
for the shell and ISM, respectively. In internal collisions, the assumption is 
that there is a fast shell and a slow shell. In the strongly shocked internal 
collision region, there may be head-on and tail collisions. In the dynam¬ 
ics of Fermi acceleration, head-on collisions are much more likely [16]. If 
Tj-ei 3> 1, then this simplifies to collisions with a stationary target. Thus, 
fiKKiE) = npj{E)a{E)np^sc/Trei, where Upj and rip ,, are the particle densi¬ 
ties for the fast and slow shell, respectively. 
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3. Results 


Our calculation gives the number of events as a function of energy, and 
a total cross section. As the collision energy rises, the KK excitation cross 
section increases and diverges further from its SM counterpart due to the 
larger influence from higher order processes EDI. The cross section for pro¬ 
ducing quark antiquark pairs is approximately an order of magnitude greater 
than for producing leptons. All cross sections are within a 95% conhdence 
interval. These Monte Carlo generated cross sections are shown in hgure 



Figure 1: PythiaS generated cross sections for KK excitations. All values are within a 
95% confidence interval. 


We calculated production rates when a shell smashes into the ISM, a 
dense molecular cloud, and a slower shell. For the shells, we used a mass 
of 10“^ Mq, which is consistent with the mass estimates in section 2.2 We 
assumed that Fermi-type acceleration generated a power-law spectrum for 
the differential particle densities in the shocks with energies from 100 to 
10® GeV and a Lorentz factor of 300. For the internal collisions, we report 
the reaction rates between two shells with the same mass, a fast shell with 




a Lorentz factor of 300, and a slower target shell with a Lorentz factor of 1. 
The KK excitation rates are reported in hgure The total reaction rates 
in cm“^s“^ were 53, 5300, and 5.3x10^^ for ISM case, dense cloud case, and 
internal collisions. 

4. Conclusions and Observational Implications 

Our results indicate that KK excitations could be signihcant as a GRB 
hreball shell collides with another medium, especially since KK mediated 
processes enjoy vastly boosted cross-sections over their SM equivalents over 
the relevant energy range. In particular, we found that collisions with the 
ISM and a dense cloud would produce KK excitations, but at a much slower 
rate than internal collisions. Despite the increase in the cross-section with 
particle energy, the excitation rate decreases as there are fewer protons at 
higher energies. 

Reference |2] found a characteristic peak and valley in the KK excitation 
rate as a function of beam energy. The suppression and boosting of the cross 
section that leads to these features is consistent with our simulations. The 
LKP should be identihable by the characteristic suppression of the number 
of events near its resonant rest mass energy (1.59TeV in our work). This 
would in turn lead to a decrease in gamma ray output from a GRB near 
and above this energy scale. Such a dark matter production signal does 
not conflict with any observational constraints, and would account for the 
signals gamma-ray excess from the Fermi and Pamela experiments. In the 
near future, improvements to the current generation of TeV instruments, 
the newly operational High-Altitude Water Gherenkov Observatory (HAWG) 
experiment, and the planned Gherenkov Telescope Array (GTA) experiment 
aim to observe GRBs in the TeV regime. In particular, GTA is expected 
to revolutionize the study of cosmic particle accelerators, and may provide 
further evidence of DM production. 

We assumed that GRB hreballs are baryon dominated, however other the¬ 
ories assume that the outflow is composed of leptons, primarily electrons. The 
resulting spectrum is then from one of or a combination of synchrotron, in¬ 
verse Gompton, and synchrotron self-Gompton. The benehts of these leptonic 
models include their relative simplicity, and the numerous adjustable param¬ 
eters for htting purposes. On the other hand, the complexity of hadronic 
models has limited their use in simulations as the codes must keep track of 
all the particles from particle showers and determine their contribution to 
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the spectrum. Despite these difficulties, there has been some recent work 
attempting to explain observations with hadronic flows [IHIEI]. Using sim¬ 
ilar methods, we hope to compute the predicted GRB spectrum taking KK 
excitations into account. 

In the future, we also intend to explore the possibility that SNR may 
be significant KK production environments. As a source of galactic cosmic 
rays, SNRs also create shocks and accelerate protons. SNRs live longer than 
GRBs and have been observed in the TeV range [32]. Furthermore, hadronic, 
photopion models have already been shown to fit the gamma-ray spectrum 
from SNRs [3^ . 
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Figure 2: Energy integrated production optical depth for KK excitations for various pa¬ 
rameters. The optical depth indicates the likelihood that a particle will undergo collisions 
resulting KK excitation while traveling across the width of a shell. The plot is a function 
of the distance between the shell and the central engine. In the case of a purely SM regime, 
this is equivalent to the production optical depth SM bosons 7 or The shell density 
is 10 ^^ cm“^ with a energy distribution that follows a power law with index a = — 2 . 
The width of the shell is 10^ cm. The Lorentz factor is 300. In the upper left panel, we 
assume multiple shells expanding with partij^ density no on the order of 10^^ cm“^. The 
upper right panel assume that the shell has slowed and the Lorentz factor is 30. Finally, 
the lower panel assume that the burst outputs a single shell with 10^^ cm“^ and Lorentz 
factor is 30. 
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Figure 3: Differential reaction rate for KK excitations. The densities are 1 cni“^ and 1000 
cm“^ for the ISM and dense molecular cloud, respectively. The isotropic energy output 
for the baryons associated with each shell is 10^^ ergs and the initial radius is 10^ cm. 
The incident particle density no was on the order of 10^^ cm“^. For internal collisions, we 
assume that the slow shell has slowed such that F^e; ~ Fy. 
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